Persistent currents in a Moebius ladder: A test of interchain coherence of interacting 
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Persistent currents in a Moebius ladder are shown to be very sensitive to the effects of intrachain 
interactions on the hopping of electrons between chains. Their periodicity as a function of flux is 
doubled for strong enough repulsive interactions because electrons cannot hop coherently between 
the chains and have to travel along the full edge of the Moebius ladder, thus encircling the flux 
twice. Mesoscopic devices that should enable one to observe these effects are proposed. 
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The problem of transport in very anisotropic systems 
has become one of the central issues in the field of 
strongly-correlated systems. Roughly speaking, it can 
be stated as follows: Consider an electronic system with 
transfer integrals much smaller in one direction than in 
the other(s). Is it possible, and under which conditions, 
that transport is coherent in the highly conducting direc- 
tion^) and incoherent in the other direction at tempera- 
tures much smaller than the smallest transfer integrals? 
This question has recently been raised in several con- 
texts, the most prominent examples being the high-T c 
superconductors for quasi-2D systems and the or- 

ganic superconductors, the Bechgaard salts, for quasi-lD 
systems |3|,[|]. This property is clearly inconsistent with 
Fermi Liquid theory, and the best candidates to describe 
such a behaviour are models of strongly-correlated elec- 
trons. Metallic behaviour is characterized by the pres- 
ence of a Drude peak in the optical conductivity at zero 
temperature, and what one is looking for is a model for 
which there is no Drude peak in the conductivity in one 
direction, while there is one in the other (s). The optical 
conductivity is not a simple object to calculate though, 
and a number of authors concentrate on the Green's func- 
tions and look for models with no pole in the Green's 
function in one direction and one in the other(s). Both 
properties are very likely related, although this is not 
totally clear due to vertex corrections. 

So far, most results have been obtained for quasi- 
lD systems, i.e. systems of coupled Luttinger liquids. 
Renormalization Group (RG) arguments || suggest that 
there will be a pole in the transverse Green's function as 
long as the Luttinger liquid exponent a is smaller than 1, 
while another approach drawing an analogy to the prob- 
lem of coherence in a two-level system coupled to a bath 
suggests that the pole in the transverse Green's func- 
tions might disappear much earlier |3j . Evidence in favor 
of the latter picture based on the analysis of the angular 
dependence of the magnetoresistance of the Bechgaard 
salt (TMTSF) 2 PF 6 has been presented |{§. However, 
the issue remains controversial, and more direct evidence 



of the effect of correlations on the coherence of motion 
perpendicular to the chains would be highly desirable. In 
the case of two chains, coherence shows up as a splitting 
of the bonding and antibonding bands, which could in 
principle be measured in angular resolved photoemission 
experiments. However, considering the difficulty in inter- 
preting photoemission experiments on low dimensional 
systems, and the need to vary parameters like a or the 
perpendicular hopping, this prospect remains remote to 
say the least. 




FIG. 1. Schematic diagram of a Moebius ladder pierced by 
a magnetic flux (7ic/e)$. 

In this paper, we propose a new approach to study the 
problem of interchain coherence based on the analysis of 
persistent currents in a Moebius ladder (see Fig. 1). If 
electrons can hop coherently from one edge to the other, 
the groundstate energy and the persistent current will 
be periodic functions of the flux with the usual period 
4>o = hc/e. However, if the perpendicular hopping in- 
tegral t± is switched off, electrons will have to go twice 
around the flux to reach the same site, and the period 
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will be 4>o/2. Now, for an interacting system with tj_ =/= 0, 
the period is expected to be 4>o &s long as coherent mo- 
tion between the chains is possible, and to become <j)o/2 
when interactions are strong enough to prevent any co- 
herent motion between the chains. Thus, measuring the 
periodicity of the persistent current should provide direct 
information on the coherence of interchain hopping. 
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FIG. 2. AE for a non-interacting, quarter filled system: a) 
as a function of t±; b) as a function of the number of rungs 
N. 

This picture turns out to be essentially true, although 
the analysis is not as straightforward as one might hope 
due to rather interesting finite-size effects. Let us start 
with a careful analysis of the non-interacting case. For 
simplicity, we consider spinless fermions throughout this 
paper. The spin is not expected to play any significant 
role for this property, a point of view confirmed by pre- 
liminary results we have obtained for fermions with spin. 
Spinless fermions on a Moebius ladder with N rungs 
pierced by a magnetic flux (7ic/e)$ can be described by 
a one-dimensional periodic Hamiltonian 

2N 



H=-tJ2^ /N 4c> +1 + h.c.) 
i=\ 

j 2N 
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with the usual convention c i+2 N = c i- This Hamiltonian 
is readily diagonalizcd by a Fourier transform, and the 
dispersion reads 



e k = -2tcos(ka + $/N) - t± cos(Nka) 



(2) 



with k = p(2ir/2Na), p integer. The essential ladder 
structure of the Moebius ladder is contained in this ex- 
pression because cos(Nka) = cos(pn) = +1 if p is even 



and — 1 if p is odd. The system thus consists of bond- 
ing and antibonding bands with the usual dispersions 
e k = ~2tcos(ka + &/N) ± t±. The difference with a 
standard ladder is that the wavevectors k = p(2n/2Na) 
are restricted to even and odd values of p for the bonding 
and antibonding bands respectively. 

Let us consider the periodicity of the groundstate en- 
ergy of such a system as a function of If t± is large 
enough, all the fermions are in the bonding band, and the 
total energy reads E($) = Y, P h 2t cos(pir/N + $/N) - 
t±], where the sum over p is restricted to even integers 
chosen to give the lowest energy for a given <!>. This func- 
tion is clearly periodic in <I> with period 2n. If t± = 0, the 
bonding and antibonding bands form a single band, and 
the total energy reads E($) = -2t ^ cos(pn/N + $/N), 
where the integers p can now be both even and odd, and 
the periodicity of this function is it. For intermediate val- 
ues of t± , the periodicity is 2n except for specific values of 
t± where it is it. The number of such points scales with 
the number of particles in the system. To understand 
this, let us consider the quantity 



AE = (-i)^N[E(<f> = tt) - E(& = 0)] 



(3) 



where n is the number of particles. In the following, n is 
restricted to odd values to have a non degenerate ground- 
state which makes the analysis slightly simpler, although 
the results are essentially equivalent for an even number 
of particles, and the factor (— 1)~3~~ has beeen included 
to insure that AE is always positive in the limit t± — > 0. 
Besides, as for the curvature at $ = 0, which gives the 
Drude weight, one has to multiply by N for ID systems 
to get non- vanishing results in the thermodynamic limit. 
A typical example of the behaviour of AE with t± is de- 
picted in Fig. 2a. It is clear from the dispersion of Eq. 
(2) that both £($ = tt) and £($ = 0), and hence AE, 
are piece- wise linear functions of t±. The slope of AE 
changes each time a pair of particles goes from the anti- 
bonding band to the bonding band, which we know has 
to occur because the antibonding band is certainly empty 
for large enough t±. Now these transitions do not occur 
for the same values of t± for $ = n and $ = 0. As a 
result, the slope of AE alternates between 2N and — 2N 
until the antibonding band is empty, in which case it is of 
course 0. Between changes of slopes, AE vanishes once, 
which corresponds to points where the periodicity is n, 
like for t± = 0. Since the particles change bands in pairs, 
the number of such points is essentially half the num- 
ber of particles (more precisely 2n + 1 including t± = 
for 4n + 1 particles, and similar formulae for other fill- 
ings). The same effect shows up in the behaviour of AE 
with N: For intermediate values of t±, AE oscillates as 
a function of size with extrema each time the difference 
in the particle number between the bonding and the an- 
tibonding bands increases by 2 (see Fig. 2b) . The points 
lie on a piece-wise linear curve with a slope alternating 
between 2t± and — 2t±. As a consequence, the periodic- 
ity of the groundstate energy for a given t± and for non 
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interacting fermions is not a well denned quantity in the 
thermodynamic limit. 




FIG. 3. AE as a function of tj_ for V = 2.5 (a = 0.36) and 
V — 5 (a — 1) for a quarter-filled system with 18 rungs. All 
energies are in units of t. 



duced by interactions and changes sign between V/t = 5 
and V/t — 6, which means that AE decreases with N 
(see Fig. 4a). Given the absence of oscillations in AE 
as a function of t± for these values of the interaction, it 
is natural to assume that AE does not change sign as a 
function of N either. This leads us to the conclusion that 
AE goes to in the thermodynamic limit when a > 1 
||. Significant effects are already present for a < 1, 
however. Let us consider for instance t± = 0.4, for which 
AE changes sign between N — 10 and N — 14 in the 
non-interacting case (see Fig. 4b). For the sizes we can 
reach, AE already does not change sign for V/t = 2, and 
the curve is very flat for V/t as small as 3 (a — 0.46). 
It is in fact possible that even in that case AE goes to 
with damped oscillations. Obtaining numerical results 
on larger systems by other methods, e.g. Densitiy Matrix 
Renormalization Group, would be very useful to check 
this point. Let us note that oscillations would also be 
present for a regular ladder. However the periodicity 
remains equal to 0o even if t± is switched off. So the 
Moebius geometry is essential to observe the effects of 
interactions. 



Let us now consider the effect of intrachain interac- 
tions on the periodicity of the groundstate energy To 
describe systems with large values of the Luttinger liq- 
uid exponent a, we consider an interaction term of the 
form 

/ 2V V \ 

H^t = y J [ Vnin i+1 + —niUi+2 + —nin i+3 j (4) 

For a quarter-filled system, the exponent a of this model 
has already been calculated with standard techniques , 
and it reaches the value 1 for V/t — 5. Using Lanczos 
technique, we have calculated the dependence of AE on 
t± for different values of V/t and different sizes up to 
36 sites (N = 18). Typical results are shown in Fig. 3. 
The effects of intrachain repulsion are rather dramatic. 
The first oscillation quickly becomes the dominant one, 
and the first value of t± where AE vanishes increases 
significantly: The curve is already strongly affected for 
V/t = 2.5 (Fig. 3) with respect to the non-interacting 
case (Fig. 2a). But more importantly, the oscillations 
disappear altogether for V/t — 5, i.e. a = 1. This crit- 
ical value turns out to be independent of the size. The 
origin of the oscillations being that particles go from the 
antibonding band to the bonding band, this means that 
these concepts have lost their meaning when a is big 
enough. In other words, t± is no longer able to produce 
two separate bands in the low energy spectrum of the 
system. 

Our discussion of the dependence of AE on N for in- 
teracting systems is limited by the maximum size we 
can handle with Lanczos, namely 36 sites (N = 18) at 
quarter-filling. For t± = 0.1, we are limited to the first 
linear section of Fig. 2b where the slope is equal to 2t± 
in the non-interacting case. The slope is considerably re- 
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FIG. 4. AE as a function of the number of rungs for 
different values of V and two values of t±: a) t± = 0.1; b) 
t± = 0.4. All energies are in units of t. 



To make contact with current theories of the effect 
of interactions on interchain hopping, let us first note 
that the present results strongly suggest that the system 
has lost any memory of the splitting between bonding 
and antibonding bands when a = 1. This is reminiscent 
of the RG result that t± is an irrelevant perturbation 
when a > 1, but it is much stronger: RG arguments 
are limited to infinitesimal values of t ± , while our results 
show that large values of t± are still unable to produce 
a difference between bonding and antibonding states if 
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a > 1. In fact, we believe that the results of Fig. 3 
are the first direct numerical evidence of a strong effect 
of interactions on hopping between chains because they 
exhibit a qualitative change without having to go to the 
thermodynamic limit. Besides, our results concerning 
the dependence of AE on N show very dramatic effects 
for relatively small interactions. These conclusions agree 
qualitatively with those of Refs. Q and ||, and they are 
consistent with previous numerical studies of the spectral 
functions 



h = (-!)' 



eAE 
2HN 



(5) 




FIG. 5. Possible realization of a Moebius ring using a wide 
quantum wire. At low densities, electrons will form Luttinger 
liquids along the edges of the wire. Interedge tunneling is 
controlled by gate gl, while the crossover between the two 
edges is controlled by gate g2. 

TThe main advantage of the present approach is that 
the quantities we calculate are not spectral functions, 
whose low-energy behaviour is very difficult to mea- 
sure with today's resolution, but thermodynamic quanti- 
ties that should be accessible in experiments on meso- 
scopic systems. A Moebius ring with controllable in- 
teredge tunneling could in principle be fabricated in a 
GaAs heterostructure: At sufficiently low carrier densi- 
ties, electrons in a wide quantum wire will form Lut- 
tinger liquids along the two edges, and interedge tun- 
neling could be controlled via an external gate gl (see 
Fig. 5). The crossover between the two edges could be 
provided by a ballistic cross, controlled by gate g2, as 
illustrated in Fig. 5. The persistent current in such a 
structure, which is related to the ground state energy by 
/($) = —(e/h)dE/d<f>, could be determined experimen- 
tally via magnetization measurements ]ic|j. can be 
expressed as a Fourier series, /(<&) = Xm=i sin(n$), 
and such experiments typically measure the first two har- 
monics, I\ and li- The first harmonic is in fact propor- 
tional to AE, 



provided the higher odd harmonics can be neglected 
(which is always the case experimentally). ^From the 
above discussion of AE, it is clear that I\ should be a 
highly sensitive function of tj_ and of the carrier density 
in the coherent tunneling regime, while I\ = in the 
absence of coherent interedge tunneling. Experiments 
on suitably fabricated mesoscopic systems, in which the 
predicted finite-size effects could be directly measured, 
are thus good candidates to test theories of the effect of 
interactions on interchain tunneling. 

In conclusion, we have shown that interactions between 
electrons have dramatic effects on their ability to hop co- 
herently from one chain to the other by studying the 
flux dependence of the ground-state energy - or equiva- 
lently of the persistent currents - in a Moebius ladder. 
The absence of oscillations as a function of t± for strong 
enough interactions is, we believe, the best evidence of 
an interaction-induced destruction of interchain coher- 
ence obtained so far with numerical simulations. Further 
work along these lines, either numerically by studying 
larger systems, or experimentally by measuring persis- 
tent currents in appropriate mesoscopic devices, should 
be a promising area for future research. 
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